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We report on the modification of the g-factor by an in-plane electric field in an In 0 . 031 Ga 0.969 As 
epilayer. We performed external magnetic field scans of the Kerr rotation of the InGaAs film in 
order to independently determine the g-factor and the spin-orbit fields. The g-factor increases from 
—0.4473 ± 0.0001 at 0 V/cm to —0.4419 ± 0.0001 at 25 V/cm applied along the [110] crystal axis. 
In addition, spatially-resolved spin measurements show a g-factor dependence on diffusive velocity. 
The change in g-factor with electric field can have a large effect on the determination of the internal 
spin-orbit and nuclear fields from Larmor precession frequency measurements. 


The field of spintronics is based on the manipulation of 
spins by electric and magnetic fields. The strength of the 
interaction of a spin with a magnetic field is determined 
by the g-factor. Therefore, local manipulation of the g- 
factor allows for local control of spins. Electron spin reso¬ 
nance (ESR) was shown to induce spin flips in a quantum 
dot when the frequency of an oscillating magnetic field 
matched the resonance frequency of the spin in a per¬ 
pendicular static magnetic field, g^s-Bstat = hf OB <A As 
this technique requires high magnitude, high frequency 
magnetic fields, alternative methods were developed that 
instead used an alternating electric field to produce the 
resonance, either by creating a oscillating spin-orbit (SO) 
field^ or by electrically modulating the g-tensor, called g- 
tensor modulation resonance (g-TMR)^. The use of an 
electric field also allows for more local manipulation than 
that possible with magnetic fields^. 

g-TMR requires electrical control of the g-factor, which 
has been demonstrated in quantum wells (QWs)£ and 
quantum dots (QDs)£ but not in bulk materials. In QWs 
and QDs, the change in the g-factor can be attributed to 
the shifting of the wave function into the barrier, which 
has a different g-factor, by an electric field. In this paper, 
we demonstrate electrical control of the g-factor in a bulk 
In 0 .031 Gao. 969 As epilayer. We characterize the g-factor 
dependence on the in-plane electric field and drift velocity 
using magnetic field and spatially-resolved pump-probe 
Kerr rotation spectroscopy. Our results for bulk InGaAs 
demonstrate that the electric field dependence here must 
be due to a different mechanism than the wave function 
shift and distinguish changes in the g-factor from changes 
in the spin-orbit field. A change in the g-factor with in¬ 
plane electric field was recently reported in a QW—. 

Understanding the change in the g-factor as a func¬ 
tion of the electric field is also important to correctly 
determine the magnitude of the internal effective mag¬ 
netic fields produced by spin-orbit coupling and nuclear 
spin polarization. Fits for the internal fields from time- 
resolved Kerr rotation dataA— assume that the g-factor 
is constant in order to separate it from the internal fields. 
In Ref. [ll|, the g-factor is known not to change, but the 
precession frequency changes due to a change in the in¬ 
ternal fields by an electric field. However, this is referred 


to as an effective tuning of the g-factor, using the rela¬ 
tionship /7eff(B)B ex t = fj (B e xt + Bj nt (B)). 

Measurements were performed on a 500 nm thick 
Si-doped Ino.o 3 iGao. 96 gAs epilayer grown by molecular 
beam epitaxy on a (001) GaAs substrate. At room tem¬ 
perature, the carrier density was measured to be 1.45 
x 10 16 cm -3 and the mobility was 4870 cm 2 /Vs. Two 
channels 400 /im long and 100 //in wide were patterned 
and ohmic contacts were deposited to apply an in-plane 
voltage either along the [110] and [110] directions. 

The g-factor of the spins is measured using a Kerr rota¬ 
tion setup. A mode-locked TitaniunnSapphire laser, with 
a 76 MHz repetition rate, tuned to 839.43 nm, is split 
into pump and probe beams. A mechanical delay line is 
used to control the temporal separation of the pump and 
probe. For our magnetic field-dependent measurements, 
the temporal separation of the pump and probe is set to 
13 ns. In order to induce a spin polarization in the sam¬ 
ple according to the optical selection rules^, the pump 
is circularly polarized. The Kerr rotation of the linearly 
polarized probe is measured using a balanced photodiode 
bridge. The pump and probe are modulated by a pho¬ 
toelastic modulator and an optical chopper respectively 
for cascaded lock-in detection. The spatial separation of 
the pump and probe is controlled with a scanning mirror. 

For photoluminescence (PL) measurements, the sam¬ 
ple is excited with a 633 nm HeNe laser. The spectrum is 
analyzed using a spectrometer with a silicon CCD. Both 
Kerr rotation and PL measurements are performed at 30 
K, unless otherwise noted. 

Both time-resolved and magnetic field-dependent 
Kerr rotation measurements can be described by the 
equation^: 

8 k (At,B ext ) = ^ ^ e -(At+n tr ep)/T 2 * cos [^£ Bext(A / +n / rep )] 
n 

, (1) 

where A is the Kerr rotation amplitude, At is the time 
delay of the pump and probe, t rep = 13.16 ns is the 
laser repetition rate, T 2 * is the spin lifetime, g is the elec¬ 
tron g-factor, /is is the Bohr magneton, Ti is the reduced 
Planck’s constant, B ex t is the magnitude of the external 
magnetic field, and n is the pulse number. 
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When an in-plane voltage is applied across the sample, 
the spin packets created by subsequent pump pulses will 
be spatially separated due to carrier drift. Equation |T| is 
modified to account for a spatially- and time-dependent 
amplitude A n (x,At). Furthermore, the spins experience 
a momentum-dependent spin-orbit splitting that acts like 
an effective internal magnetic field, requiring that the 
external magnetic field in Eq. |T] be replaced by a vector 
sum of the internal and external magnetic fields^: 


@k (B ex t , x) ^ ' A n (x, At) x 

n 

- \j ^-Bext + -B||^ + B 2 ^_(At + nt rep ) 


99B 


( 2 ) 


where B\\ and B± are the components of the internal 
field parallel and perpendicular to the external field, re¬ 
spectively. 

Varying At at a fixed B ext results in oscillations that 
are periodic in J (^B ext + B\\ j + B\. If instead At 


is fixed and B ext is varied, then the oscillations are peri¬ 


odic in tLtLS-A t. 

n 

The components of the spin-orbit field parallel and 
perpendicular to the external field shift the curve and 
decrease the magnitude of the center peak respectively 
(Fig. HJi,b), whereas the g-factor changes the frequency of 
the peaks (Fig. [l}c)- In this way, unlike in time-resolved 
measurements, we can separately determine the g-factor 
and the spin-orbit field. 

Magnetic field scans were performed at different pump- 
probe separations, with the external magnetic field par¬ 
allel to the applied voltage. The amplitude as a function 
of this separation (Fig. [2jt,b,c) shows the profile of the 
spin packet. The difference in amplitude between the 
+2 V and the -2 V scans is due to drift of various fac¬ 
tors over the course of the measurements, including laser 
power, vertical alignment of the pump, and the balance of 
the photodiode bridge. However, the absolute amplitude 
doesn’t affect the subsequent fits. 

The average drift velocity ( Vd ) of the electrons can be 
calculated as a function of voltage from the center po¬ 
sition ( x c ) of the amplitude and the known temporal 
separation of the pump and probe pulses. From mea¬ 
surements of the drift velocity as a function of applied 
electric field, we determine the effective mobility, using 
the equation /x e ff = i’d/E ap p ii ec j, to be 5110±70 cm 2 /(V s) 
along the [110] direction and 5800 ±200 cm 2 /(V s) along 
the [110] direction. The difference in mobilities along 
the two directions comes from a discrepancy between the 
total applied voltage and the actual voltage across the 
channel, as there is likely also a voltage drop across the 
contacts that is different for different contacts. As such, 
all electric field dependent measurements are presented 
in terms of the drift velocity, which we can measure di¬ 
rectly. 
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FIG. 1. (a) Kerr rotation data as a function of the external 
magnetic field for applied voltages of 0 V (black), 1 V (red), 
and 2 V (blue) V fit to Eq. 2. The fit values were then plugged 
back into Eq. 2 and plotted varying only the spin-orbit field 
(b) and only the g-factor (c) as a function of voltage to show 
the separate dependences on these two variables. 


SO field strength measurements were performed as in 
Ref. [lU (Fig. 2d,e,f). The internal fields were found to 
be perpendicular to the direction of the current for both 
crystal axes, as is expected for these directions. The 
spin-orbit field proportionality constants were found to 
be 1.96 ± 0.02 mT ns pm -1 along [110] and 1.02 ± 0.02 
mT ns /xm _1 along [110] (Fig. 3b). 

The g-factor was extracted from the external magnetic 
field scan at each pump-probe spatial separation (Fig. 
[2k,h,i). The position dependence of the g-factor is likely 
related to the velocity-dependent spatial separation of 
the spins, indicating that there is some diffusive velocity 
dependence to the g-factor. 

Since the probe beam has some width w p , the mea¬ 
sured g-factor at a given pump-probe separation is an 
average of the g-factor of the electrons within the probe 
spot. Given the probe and spin packet widths, we can 
deconvolute the effects of the finite probe size from the 
measured g-factor to get the corrected g-factor. Since 
the probe beam and spin packet both have a Gaussian 
profile, we can write the measured g-factor as a function 
of the corrected g-factor—: 
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(s-ttrf.) 2 (x-xp ) 2 

rf° g(a;)e 2 ”^ e 2 ”r da: 

5meas(a:p) = (a ._ Xrf) 2 (x _ aip )2 ( 3 ) 

/ e a e p da; 

where Wd {w p ) is the width of the spin packet (probe) 
and Xd (x p ) is the position of the spin packet (probe). 
Since the measured g-factor has a linear dependence on 
position, we can assume that the corrected g-factor has 
the same dependence. Using this assumption, we can 
integrate Eq. [3] and rearrange to solve for the corrected 
g-factor: 


g(4) — 9meas{x') “f“ ( J (tineas(**0 f?meas(*£rf)) (d) 

\W d J 

Both the measured g-factor and the corrected g-factor, 
calculated in this way, are shown in Fig. 2g,h,i, using 
w p = 17.5 ±0.5 fin i and Wd determined from a decon¬ 
volution of the probe with the amplitude profiles in Fig. 
la,b,c. One can see that the deconvolution of the g-factor 
increases the change in the g-factor as a function of po¬ 
sition. Since the pump is at x = 0 fi m, spins that are 
farther from the pump location have a higher diffusive 
velocity. Thus, we can conclude that faster moving elec¬ 
trons have a g-factor that is larger in magnitude. 

The g-factor at the center of the spin packet is plotted 
versus the drift velocity Vd (Fig. [3k). Parabolic fits are 
shown as a guide to the eye. However, a numerical deriva¬ 
tive shows that this fit is less accurate for larger drift 
velocities. Measurements conducted along the [110] and 
[110] crystal axes show similar curvatures: (1.04 ± 0.7) 
xl0- 3 /(^f) 2 for [110] and (1.23 ± 0.2) xl0- 3 /(^f) 2 
for [110]. On the other hand, the spin-orbit field pro¬ 
portionality constant is about twice as large along the 
[110] direction. The offset is likely due to the anisotropic 
g-factor tensor—. 

It is interesting to note that we see two opposite effects 
on the g-factor due to diffusive velocity and the applied 
electric field. Figures [2k, i show that spins that have ex¬ 
perienced a larger net displacement have a higher mag¬ 
nitude g-factor, whereas Fig. [3k shows that electrons in 
a larger electric field have a smaller magnitude g-factor. 

Ref. [13 indicates that the energy dependence of the g- 
factor in GaAs is given by g*{E) = g* + 6.3 eV” 1 • E. If 
we assume a parabolic dispersion relation for our values 
of k , then our measured dependence on the energy, taken 
from electric-field dependent measurements, is more than 
500 times larger. Furthermore, we see the opposite trend 
for electrons of greater speed. 

Measurements for both positive and negative Vd shows 
that the change is g-factor is symmetric about zero drift 
velocity (Fig. |4ja) , and thus only depends on the mag¬ 
nitude of Vd, not the direction. Temperature-dependent 
measurements of the g-factor (Fig. |3 ]d) were performed 
on the same channel in order to compare the effects of 
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FIG. 2. (a), (b), (c) Amplitude of the Kerr rotation, (d), 
(e), (f) fits of the perpendicular component of the internal 
field, and (g), (h), (i) fits of the g-factor as a function of 
pump-probe separation for -2 V, 0 V, and +2 V. The spatial 
dependence of the g-factor that indicates a diffusive velocity 
dependence. Corrected g-factor values found using Eq. 4 are 
shown in gold. 


temperature and electric field. The change in g-factor 
due to an applied voltage of 1 V at 30 K is equivalent 
to the change due to the channel heating by 9 K. From 
power dissipation calculations, we can estimate that the 
expected change in temperature of the channel due to an 
applied voltage of 1 V is 0(1O -4 ) K. Therefore, the g- 
factor dependence on voltage and temperature are likely 
distinct phenomena. 

To check whether the applied voltage was causing ex¬ 
cessive channel heating, we performed temperature and 
voltage dependent PL measurements. From 30 K to 40 
K, there was a clear shift in the energy of the peak due 
to the change in the bandgap with lattice temperature. 
However, the PL for 0 V and 2 V had no discernible shift. 
Therefore, it is unlikely that the applied voltage is caus¬ 
ing sufficient channel heating to account for the change 
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FIG. 3. (a) g-factor as a function of electric field (shown in 
terms of the drift velocity) for the [110] and [1-10] directions. 
The parabolic fits to the g-factor are a guide to the eye and 
are less accurate a higher drift velocities. The offset is likely 
due to g-factor anisotropy, (b) The internal field magnitude 
as a function of drift velocity. The slopes of the lines give the 
spin-orbit field proportionality constant. 

in the g-factor. 

We expect the electron temperature to be signifi¬ 
cantly above the lattice temperature for our range of 
applied electric fields, and one possible explanation for 
the change in g-factor is that the g-factor is dependent 
on the electron temperature. It was reported that in n- 
GaAs the electron temperature increased from 4.2 K at 0 
V/cm to 38 K at 20 V/cm and 75 K at 50 V/cm!£. The 
electron temperature was estimated from the high energy 
tail of the PL. We saw no change in the PL of our sample 
as a function of electric field. However, we expect the PL 
of our sample to have greater broadening due to it being 
an alloy. The broadening due to alloy fluctuations could 
be overwhelming any changes due to the electric field. 

The electron temperature can be estimated if the en¬ 
ergy loss rate per electron is known, based on mea¬ 
surements done on AlGaAs/GaAs heterostructures using 


Electric field (V/cm) 
-18.5 0.0 18.5 
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FIG. 4. (a) g-factor as a function of drift velocity along [1- 
10] at 30K, with a parabolic fit as a guide to the eye. (b) 
Temperature dependence of the g-factor for vd = 0. The 
change in g-factor due to an applied voltage of 1 V at 30 K 
corresponds to the change to the channel heating by 9 K. 


PLi2, Shubnikov-de Haas Oscillations^, and far-infrared 
spectroscopy^!. For an applied voltage of 2 V, we found 
the energy loss rate per electron in our samples to be 
2 x 1(T 12 W, which corresponds to an electron tempera¬ 
ture of about 50 K for a lattice temperature of 4.2 K. 

Another possible explanation is that the applied elec¬ 
tric held is modifying the wave-function of donor-bound 
electrons, causing their wave functions to spread farther 
into the semiconductor, thus changing the g-factor—. 
Calculations for a Si dopant in GaAs show that the rel¬ 
ative magnitude change in the g-factor is comparable to 
what we measure here. However, the sign of the change 
is opposite. Furthermore, for the doping density of the 
sample and at the temperatures considered here, con¬ 
tributions from donor-bound electrons are expected to 
be minimal. Therefore we can rule out modification of 
donor-bound electron wave functions as the cause of the 
g-factor modification. 

The orbital contribution of spin-orbit induced circulat¬ 
ing currents was shown to be significant in calculations of 
the g-factor in quantum dots2!. Similarly, the net drift 
velocity of the spins in our sample could be modifying 
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the orbital contribution to the g-factor. 

We have performed electric-field dependent measure¬ 
ments of the g-factor in a bulk Ino .031 Gao. 969 As epilayer 
in a manner that distinguishes between changes in the g- 
factor from changes in the spin-orbit field. Separate de¬ 
termination of these two quantities is important as their 
percent change with voltage is comparable. For exam¬ 
ple, for measurements along the [110] direction (Fig. 3), 
we found the g-factor to be —0.45268 ± 0.0003 and the 
SO-field to be 4.50 ± 0.08 mT for the largest drift ve¬ 
locity. If instead time-resolved measurements had been 
done to determine the SO field from the Larmor preces¬ 
sion frequency with an applied magnetic field of 0.2 T 
assuming the zero field value of the g-factor, the SO field 


would have been calculated to be only 2.33 ± 0.20 mT. 
However, more work is needed to develop a quantitative 
model for this phenomenon. 
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